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Abstract 

We report the first ab initio density-functional study 
of (111) screw dislocations cores in the bcc transition 
metals Mo and Ta. Our results suggest a new pic- 
ture of bcc plasticity with symmetric and compact 
dislocation cores, contrary to the presently accepted 
picture based on continuum and interatomic poten- 
tials. Core energy scales in this new picture are in 
much better agreement with the Peierls energy barri- 
ers to dislocation motion suggested by experiments. 

The microscopic origins of plasticity are far more 
complex and less well understood in bcc metals than 
in their fee and hep counterparts. For example, slip 
planes in fee and hep metals are almost invariably 
close-packed, whereas in bcc materials many slip sys- 
tems can be active. Moreover, bcc metals violate 
the Schmid law that the resistance to plastic flow is 
constant and independent of slip system and applied 
stress 0. 

Detailed, microscopic observations have estab- 
lished that in bcc metals at low temperatures, long, 
low-mobility (111) screw dislocations control the 
plasticity^], ||. Over the last four decades, the 
dominant microscopic picture of bcc plasticity in- 
volves a complex core structure for these disloca- 
tions. The key ingredient of this intricate picture 
is an extended, non-planar sessile core which must 
contract before it moves. The first such proposed 
structure respected the symmetry of the underlying 
lattice and extended over many lattice constants j|. 
More recent and currently accepted theories, based 
on interatomic potentials, predict extension over sev- 
eral lattice constants and spontaneously broken lat- 
tice symmetry 0, SI EL @|. While these models can 
explain the overall non-Schmid behavior, their pre- 
dicted magnitude for the critical stress required to 
move dislocations (Peierls stress) is uniformly too 
large by a factor of about three when compared 



to experimental yield stresses extrapolated to zero- 
temperature B [tJ. 

We take the first ab initio look at dislocation core 
structure in bcc transition metals. Although we 
study two metals with quite different mechanical be- 
havior, molybdenum and tantalum, a consistent pat- 
tern emerges from our results which, should it with- 
stand the test of time, will require rethinking the 
presently accepted picture. Specifically, we find screw 
dislocation cores with compact structures, without 
broken symmetry, and with energy scales which ap- 
pear to be in much better accord with experimental 
Peierls barriers. 

Ab initio methodology - Our ab initio calcula- 
tions for Mo and Ta are carried out within the total- 
energy plane- wave density functional pseudopotential 
approach^, using the Perdew-ZungerQ parameteri- 
zation of the Ceper ly- Alder j[o) exchange-correlation 
energy. Non-local pseudopotentials of the Kleinman- 
Bylander form fllfl are used with s, p, and d channels. 
The Mo potential is optimized according toE3] and 
the Ta potential is from||l3|. We use plane wave ba- 
sis sets with energy cutoffs of 45 Ryd for Mo and 
40 Ryd for Ta to expand the wave functions of the 
valence (outermost s and d) electrons. Calculations 
in bulk show these cutoffs to give total system en- 
ergies to within 0.01 eV/atom. We carry out elec- 
tronic minimizations using the analytically continued 
approach [fl4"| within the DFT++ formalism [p~5| . 

To gauge the reliability of the pseudopotentials, 
Table |lj displays our ab initio results for the ma- 
terials' lattice constants and those elastic moduli 
most relevant for the study of (111) screw disloca- 
tions. The tabulated moduli describe the long-range 
elastic fields of the dislocations (A"), the coupling 
of displacement gradients along the dislocation axis 
z to core-size changes in the orthogonal x,y plane 
(c X x,zz = (en + 5ci2 — 2c44)/6), and the coupling of 
core-size changes to themselves in the plane (c XXiXX = 
(cn+ci2 + 2c 4 4)/2 and (c u + 2ci2-2c 4 4)/3). 
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These results indicate that our predicted core energy 
differences should be reliable to within better than 
~ 30%, which suffices for the purposes of our study. 

Preparation of dislocation cells - The cell we 
use for dislocation studies has lattice vectors Si = 
5o[l, — 1,0], a% = 3o[l, 1, — 2], and a 3 = a[l, 1, l]/2, 
where a is the lattice constant. We call this ninety- 
atom cell the "5x3" cell in reference to the lengths 
of Si and S2, and the Burgers vectors of all of the 
dislocations in our work are along 03. Eight ^-points 
ki = k 2 = j,k 3 6 ±{jq, jq, jq, jq} sample the Bril- 
louin zone in conjunction with a non-zero electronic 
temperature of fcgT = 0.1 eV, which facilitates the 
sampling of the Fermi surface. These choices give 
total energies to within 0.01 eV/atom. 

Given the relatively small cell size, we wish to 
minimize the overall strain and the effects of peri- 
odic images. We therefore follow jL6| and employ a 
quadrupolar arrangement of dislocations (a rectan- 
gular checkerboard pattern in the Si, S2 plane). This 
ensures that dislocation interactions enter only at the 
quadrupolar level and that the net force on each core 
is zero by symmetry, thereby minimizing perturba- 
tions of core structure due to the images. As was 
found in fl^j and as we explore in detail below, we 
find very limited impact of finite-size effects on the 
cores when following this approach. 

In bcc structures, screw dislocations are known 
to have two inequivalent core configurations, termed 
"easy" and "hard" |§ @, §|. These cores can be ob- 
tained from one another by reversing the Burgers 
vector of a dislocation line while holding the line 
at a fixed position. We produce cells with either 
only easy or only hard cores in this way. To cre- 
ate atomic structures for the cores, we proceed in 
three stages. First, we begin with atomic positions 
determined from isotropic elasticity theory for our 
periodic array of dislocations. Next, we relax this 
structure to the closest local energy minimum within 
the interatomic MGPT model for Mof|. faince we 
do not have an interatomic potential for Ta and ex- 
pect similar structures in Ta and Mo[[|, we create 
suitable Ta cells by scaling the optimized MGPT Mo 
structures by the ratio of the materials' lattice con- 
stants. Finally, we perform standard ab initio atomic 
relaxations on the resulting MGPT structures until 
all ionic forces in all axial directions are less than 
0.06 eV/A. 

Extraction of core energies - The energy of a long, 
straight dislocation line with Burgers vector b is E = 
E c (r c ) + Kb 3 ln(L/r c ) per b along the hnep7[, where 
L is a large-length cutoff, and K is an elastic modulus 
(see Table QJ) computable within anisotropic elasticity 



theory p8[. The core radius r c is a short-length cutoff 
inside of which the continuum description fails and 
the discrete lattice and electronic structure of the core 
become important. E c (r c ) measures the associated 
"core energy" , which, due to severe distortions in the 
core, is most reliably calculated by ab initio methods. 

The energy of our periodic cell contains both 
the energy of four dislocation cores and the energy 
stored outside the core radii in the long-range elas- 
tic fields. To separate these contributions, we start 
with the fact that two straight dislocations at a 
distance d with equal and opposite Burgers vectors 
have an anisotropic elastic energy per b given by 
E = 2E c (r c ) + 2Kb 3 \n(d/r c ). Next, by regularizing 
the infinite sum of this logarithmically divergent pair 
interaction, we find that the energy per dislocation 
per b in our cell is given by 



E = E c {r c ) + Kb 3 
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The function A(x) contains all the effects of the in- 
finite Ewald-like sums of dislocation interactions and 
has the value A = -0.598 846 386 for our cell. Sub- 
tracting the long-range elastic contribution (the sec- 
ond term of (Q)) from the total energy, we arrive at 
the core energy E c . 

To test the feasibility of this approach, we com- 
pare E c (r c ) for the MGPT potential as extracted with 
the above procedure from cells of two different sizes: 
the 5x3 cell and the corresponding 9x5 cell. (The 
MGPT is fit to reproduce experimental elastic mod- 
uli, so K is given in Table |l|.) With the choice r c = 2b, 
Table § shows that our results, even for the 5x3 cell, 
compare quite favorably with those of@, |l9), espe- 
cially given that our 5x3 and 9x5 cells contain only 
ninety and 270 atoms respectively, whereas the cited 
works used cylindrical cells with a single dislocation 
and two thousand atoms or more. Given the suitabil- 
ity of the 5x3 cell, all ab initio results reported below 
are carried out in this cell. 

Ab initio core energies - Except for the Mo 
hard core, all the core structures relax quite read- 
ily from their MGPT configurations to their equilib- 
rium ab initio structures. The Mo hard-core con- 
figuration, however, spontaneously relaxes into easy 
cores, strongly indicating that the hard core, while 
meta-stable within MGPT by only 0.02 eV/6@, is 
not stable in density functional theory. We do not 
believe that this instability is due to finite-size ef- 
fects, which appear to be quite small for the reasons 
outlined previously. 

Table |^ compares our ab initio results to available 
MGPT results for core energies in Mo and Ta. To 



2 



t • I 

..; : Im - 



1 • t 

T * ■ • M t f 

. . "'* • v 



(a) 



4-. 



<110> 



(b) 



Figure 1: DD maps as found in the MGPT model for 
(a) easy, and (b) hard Mo dislocation cores. Dotted 
lines indicate axes of C2 symmetry of the D$ symme- 
try group. 



make comparison with the MGPT, for the unstable 
Mo hard core we evaluate the ab initio core energy 
at the optimal MGPT atomic configuration (column 
AI* in Table ||). Note that, in computing hard-easy 
core energy differences, the long-range elastic con- 
tributions cancel so that these differences are much 
better converged than the absolute core energies. 

Table | shows that the MGPT hard -easy core en- 
ergy differences are much larger than the correspond- 
ing ab initio values by approximately a factor of 
three. The accuracy of the elastic moduli of Table [|, 
combined with the high transferability of the local- 
density pseudopotential approach, indicates that this 
factor of three is not an artifact of our approxima- 
tions. We believe that the reason for this discrepancy 
is that the MGPT is less transferable. Having been 
forced to fit bulk elastic moduli and thus long-range 
distortions, the MGPT may not describe the short 
wavelength distortions in the cores with high accu- 
racy. An examination of Mo phonons along [100] 
provides poignant evidence: the MGPT frequencies 
away from the zone center are too large when com- 
pared to experimental and band-theoretic values 
and translate into spring constants that are up to 
approximately three times too large. 

The magnitude of the core energy difference has im- 
portant implications for the magnitude of the Peierls 
energy barrier and Peierls stress for the motion of 
screw dislocations in Mo and Ta. In a recent Mo 



MGPT study) 19 1, the most likely path for dislocation 
motion was identified to be the (112) direction: the 
moving dislocation core changes from easy to hard 
and back to easy as it shifts along (112). The en- 
ergy barrier was found to be 0.26 eV/6, very close to 
the MGPT hard-easy energy difference itself. The 
fact that the ab initio hard-easy energy differences 



in Mo and Ta are smaller by about a factor of three 
than the respective interatomic values suggests that 
the ab initio energy landscape for the process has a 
correspondingly smaller scale. If so, the Peierls stress 
in Mo and Ta should also be correspondingly smaller 
and in much better agreement with the values sug- 
gested by experiments. 

Dislocation core structures - Figure |l] shows differ- 
ential displacement (DD) mapsQ of the core struc- 
tures we find in our ninety-atom supercell when work- 
ing with the interatomic MGPT potential for Mo. 
Our DD maps show the atomic structure projected 
onto the (111) plane. The vector between a pair 
of atomic columns is proportional to the change in 
the [111] separation of the columns due to the pres- 
ence of the dislocations. The maps show that both 
easy and hard cores have approximate 3-fold rota- 
tional (C3) point-group symmetry about the out-of- 
page [111] axis through the center of each map. The 
small deviations from this symmetry reflect the weak- 
ness of finite-size effects in our quadrupolar cell. The 
hard core has three additional 2- fold rotational (C2) 
symmetries about the three (110) axes marked in the 
maps, increasing its point-group symmetry to the di- 
hedral group D3 which is shared by the underlying 
crystal. The easy core, however, shows a strong spon- 
taneous breaking of this symmetry: its core spreads 
along only three out of the six possible (112) direc- 
tions. Our results reproduce those of|], [l9| who em- 
ployed much larger cylindrical cells with open bound- 
aries, underscoring the suitability of our cell for de- 
termining core structure. This symmetry-breaking 
core extension is that which has been theorized to 
explain the relative immobility of screw dislocations 
and violation of the Schmid law in bcc metals. 

Figure || displays DD maps of our ab initio core 
structures. Contrary to the atomistic results, we find 
that the low-energy easy cores in Mo and Ta have 
full D 3 symmetry and do not spread along the (112) 
directions. Combining this with the above results 
concerning core energetics, we have two examples for 
which our pseudopotentials are sufficiently accurate 
to disprove the conventional wisdom that generic bcc 
metallic systems require broken symmetry in the core 
to explain the observed immobility of screw disloca- 
tions. 

Turning to the hard core structures, the ab initio 
resuts for Ta show a significant distortion when com- 
pared to the atomistic core (contrast Figure [l]b and 
Figure |^b). As the ab initio Mo hard core was un- 
stable, we believe that this distortion of the Ta hard 
core suggests that this core is much less stable within 
density functional theory than in the atomistic po- 
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Figure 2: DD maps of the ab initio dislocation cores: 
(a) Ta easy, (b) Ta hard; and (c) Mo easy. Dotted 
lines indicate axes of C2 symmetry of the D3 symme- 
try group. 



tentials. 

To complete the specification of the three- 
dimensional ab initio structure of easy cores in Mo 
and Ta, Figure ^ presents maps of the atomic dis- 
placement in the (111) plane. The small atomic 
shifts, which are due entirely to anisotropic effects, 
are shown as in-plane vectors centered on the bulk 
atomic positions and magnified by a factor of fifty. To 
reduce noise in the figure, before plotting we perform 
C3 symmetrization of the atomic positions about the 
[111] axis passing through the center of the figure. 
As all the dislocation cores in our study have a mini- 
mum of C3 symmetry, this procedure does not hinder 
the identification of possible spontaneous breaking of 
the larger D3 symmetry group. Our maps indicate 
that the easy cores in both Mo and Ta have full D3 
symmetry. 

Recent high-resolution electron microscopy explo- 
rations of the symmetry of dislocations in Mo have 
focused on the small shifts in the (111) plane of 
columns of atoms along [111][^0|. This pioneering 
work reports in-plane displacements extending over 
a range much greater than the corresponding MGPT 
results and also much greater than what we find ab 
initio. In pp[| this is attributed to possible stresses 
from thickness variations and foil bending. We be- 
lieve this makes study of the internal structure of the 



Figure 3: Planar displacement maps of the ab initio 
(a) Mo easy, and (b) Ta easy cores. Vectors show 
in-plane (111) atomic shifts and have been magnified 
by a factor of fifty. Dotted lines indicate the C2 axes 
of the D3 symmetry group. 



core difficult, and that cleaner experimental results 
are required to resolve the nature of the symmetry of 
the core and its extension. 

In conclusion, our first principles results show no 
preferential spreading or symmetry breaking of the 
dislocation cores and exhibit an energy landscape 
with the proper scales to explain the observed immo- 
bility of dislocations. Atomistic models which demon- 
strate core spreading and symmetry breaking, both of 
which tend to reduce the mobility of the dislocations, 
are well-known to over-predict the Peierls stress. The 
combination of these two sets of observations argues 
strongly in favor of much more compact and sym- 
metric bec screw dislocation cores than presently be- 
lieved. 
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3.10 


3.15 


-1.6% 


3.25 


3.30 


-1.5% 
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1.60 


1.36 


18% 


0.65 


0.62 


5% 


C-xx,zz 


2.17 


1.91 


14% 


1.72 


1.39 


24% 


Cxx,xx 


5.48 


4.25 


29% 


3.02 


2.98 


1.3% 


Cxx^yy 


2.21 


1.77 


25% 


1.72 


1.49 


15% 



Table 1: Lattice constants a (A) and elastic moduli (Mbar) for Mo and Ta based on ab initio pseudopotentials 
(AI) and experiments (Expt)fl3||. 



E c (eV/6) 


5x3 


9x5 


Cylindrical|, [jj) 


hard 


2.57 


2.57 


2.66 


easy 


2.35 


2.31 


2.42 


A 


0.22 


0.26 


0.24 



Table 2: Core energies for r c = 2b as predicted by the MGPT model, easy and hard refer to different core 
configurations. A is the hard-easy core energy difference. 





Mo 


Ta 


E c (eV/6) 


MGPT 


AI* 


AI 


MGPT 


6] 


AI 


hard 


2.57 


2.94 






0.91 


easy 


2.35 


2.86 


2.64 




0.86 


A 


0.22 


0.08 




0.14 


0.05 



Table 3: Core energies for r c = 2b for fully relaxed ab initio cores (AI) and interatomic (MGPT) cores in 
the 5x3 cell. AI* refers to ab initio core energies computed based on relaxed MGPT configurations as the 
ab initio Mo hard core is unstable. Ref. Q only reported A for Ta. 
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